Making Massless Dirac Fermions from a Patterned 
Two-Dimensional Electron Gas 
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Analysis ol the electronic structure of an ordinary two-dimensional electron gas (2DEG) under an 
appropriate external periodic potential of hexagonal symmetry reveals that massless Dirac fermions 
are generated near the corners of the supercell Brillouin zone. The required potential parameters 
are found to be achievable under or close to laboratory conditions. Moreover, the group velocity is 
tunable by changing either the effective mass of the 2DEG or the lattice parameter of the external 
potential, and it is insensitive to the potential amplitude. The finding should provide a new class of 
systems other than graphene for investigating and exploiting massless Dirac fermions using 2DEGs 
in semiconductors. 



Graphenei^i 3 -^, a honeycomb lattice of carbon atoms, 
is composed of two equivalent sublattices of atoms. The 
dynamics of the low-energy charge carriers in graphene 
may be described to a high degree of accuracy by a mass- 
less Dirac equation with a two-component pseudospin ba- 
sis which denotes the amplitudes of the electronic states 
on these two sublattices. The quasiparticles have a lin- 
ear energy dispersion near the corners K and K' (the 
Dirac points) of the hexagonal Brillouin zone^^i^. Con- 
sequently, the density of states (DOS) varies linearly and 
vanishes at the Dirac point energy. The sublattice de- 
gree of freedom of the wavefunctions is given by a pseu- 
dospin vector that is either parallel or anti-parallel to the 
wavevector measured from the Dirac point, giving rise to 
a chirality being 1 or —1, respectively-^^. These two fun- 
damental properties of graphene, linear energy dispersion 
and the chiral nature of the quasiparticles, result in in- 
teresting phenomenon such as half-integer quantum Hall 
effect 2,3 , Klein paradox^, and suppression of backscat- 
tcring 6 '-ii2, as well as some novel predicted properties 
such as electron supercollimation in graphene superlat- 
ticesli^s. 

As a possible realization of another two-dimensional 
(2D) massless Dirac particle system, theoretical studies 
on the physical properties of particles in optical honey- 
comb lattices^ have been performe d 15 ' 16 i 17 i 18 i 19 . The 
behaviors of ultra-cold atoms in a honeycomb lattice po- 
tential were considered, in principle, to be equivalent to 
those of the low-energy charge carriers in graphene^. 

In this Letter, we propose a different practical scheme 
for generating massless Dirac fermions. We show through 
exact numerical calculations within an independent par- 
ticle picture that applying an appropriate nanometer- 
scale periodic potential with hexagonal symmetry onto 
conventional two-dimensional electron gases (2DEGs) 
will generate massless Dirac fermions at the corners of 
the supercell Brillouin zone (SBZ). We find that the po- 
tential configurations needed should be within or close to 
current laboratory capabilities, and this approach could 
benefit from the highly developed experimental tech- 
niques of 2DEG physics 2 - including recent advances in 
self-assembly nanostructure s 21 ' 22 i 23 . 



We moreover find that the band velocity and the en- 
ergy window within which the dispersion is linear may 
be varied by changing the superlattice parameters or the 
effective mass of the host 2DEG, thus providing a dif- 
ferent class of massless Dirac fermion systems for study 
and application. Interestingly, the amplitude of the peri- 
odic potential does not affect the band velocity about the 
Dirac points. But, if the external periodic potential is too 
weak, there is no energy window within which the DOS 
vanishes linearly. The linear energy dispersion and the 
chiral nature of states around the Dirac points of these 
2DEG superlattices are found to be identical to those of 
graphene. Also, the associated up and down pseudospin 
states naturally correspond to states localizing, respec- 
tively, on two equivalent sublattice sites formed by the 
superlattice potential. 

In order to investigate the properties of charge carriers 
in 2DEGs under an external periodic potential, we work 
within the one-electron framework with realistic material 
parameters. However, the point of this Letter is not the 
methodology itself, but the new idea and its feasibility. 
The study makes a connection between two of the largest 
and most active fields in condensed matter physics and 
nanoscience these days, graphene and 2DEGs, and pro- 
vides experimentalists working on 2DEGs with a previ- 
ously unexploited way of generating a new class of 2D 
massless Dirac fermions and of tuning their properties. 

Let us consider a 2DEG with E(p) = p 2 /2m* where 
m* is the electron effective mass. This is a good approx- 
imation to the energy dispersion of the lowest conduc- 
tion band in diamond- or zinc-blende- type semiconductor 
quantum wells, which have effective masses ranging from 
m* = 0.02 m e ~ 0.17 m e where m e is the free electron 
mass 24 . We shall consider explicitly two cases in the nu- 
merical calculations: m* = 0.02 m e and m* = 0.04 m e . 

Figure QJa) shows the muffin-tin periodic potential 
considered in our numerical calculations, whose value is 
Uq (> 0) in a triangular array of disks of diameter d and 
zero outside. Figure [IJb) is the corresponding Brillouin 
zone. The muffin-tin form is chosen for ease of discus- 
sion; the conclusions presented here are generally valid 
for any hexagonal potential. 
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FIG. 1: (a) A muffin-tin type of hexagonal periodic potential 
with a spatial period L. The potential is Uo (> 0) inside the 
gray disks with diameter d and zero outside, (b) The Brillouin 
zone of hexagonal lattice in (a) . 
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FIG. 2: (color online). Calculated energy bandstructure of 
the lowest two bands of a hexagonal 2DEG superlattice shown 
in Fig. []Ja). (a) m* — 0.02 m e and L = 20 nm, (b) m* = 
0.02 m e and L — 40 nm, (c) m* = 0.04 m e and L — 20 nm, 
and (d) m* = 0.04 m e and L = 40 nm. The diameter of the 
disks d is set to d = 0.663 L (see text). Solid red, dashed 
green, and dash-dotted blue lines show results for Uo equal 
to 200 meV, 100 meV, and 50 meV, respectively. The Dirac 
point energy (i. e. , the energy at the crossing of the two bands 
at K) is set to zero. 



We shall first discuss our numerical results and later 
consider the approximate analytic solutions. Figure [3] 
shows the calculated bandstructures of the lowest two 
bands for several hexagonal 2DEG superlattices with dif- 
ferent effective mass m* , lattice parameter L (with poten- 
tial barrier diameter d — 0.663 L^), and barrier height 
Uq. As the barrier height is decreased, the energy win- 
dow within which the energy dispersion is linear is re- 
duced (Fig. (The potential barrier heights used in 
our calculations are typical of values employed in confin- 
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FIG. 3: The DOS of a triangular 2DEG superlattice with 
m* = 0.02 m e , U = 200 meV, L = 20 nm and d = 13.3 nm. 
(The zero of energy is set at the Dirac point.) The charge 
density needed to fill the conduction band up to the Dirac 



point energy is 5.7 x 10 cm 



ing 2DEGs 20 .) However, the group velocity at the Dirac 
point is insensitive to Uo (Fig. [5]). But, as m* or L is 
increased, the group velocity decreases (Fig. [2]). 

Figure [3] shows the DOS of a hexagonal 2DEG super- 
lattice with m* = 0.02 m e , U = 200 meV, L = 20 nm 
and d = 13.3 nm. The DOS has a linear behavior around 
the Dirac point energy within a ~ 30 meV energy win- 
dow. The charge density required to dope the system to 
reach the Dirac point energy is 5.7 x 10 11 cm -2 , which 
is in the range of typical value in 2DEG studies 2 ^, and 
may be tuned by applying a gate voltage or by light illu- 
mination 2 ^. 

With the above results established from the numerical 
calculations, to gain further insight, we now present ana- 
lytical expressions for the energy dispersion relation and 
wavefunctions around the SBZ corners obtained from de- 
generate perturbation theory. We concentrate on states 
with wavevector k + K near the K point in the SBZ, 
i.e., |k| <§; |K|. Let us set the energy of the empty lat- 
tice bandstructure at the K point to zero, define W as 
the Fourier component of the periodic potential connect- 
ing 1 ->■ 2, 2 -> 3 and 3 -> 1 in Fig.[T(b) 25 , and denote 
vq as the group velocity of the electron state at the K 
point of the 2DEG before applying the periodic poten- 



tial. [For E{p) = p 2 /2r 



vo 



fiK/m*. However, the 



derivation below is not confined to a quadratic energy 
dispersion for the original 2DEG.] Due to the inversion 
symmetry of the system considered here, W is real. The 
wavefunction V'k(r) may be approximately expressed as 
a linear combination of three planewave states 



^k(r) 



[cj exp (i(Ki + k) • r) 



-c 2 exp {i{K 2 + k) • r) + c 3 exp (i(K 3 + k) • r) 



(1) 



where A c is the area of the 2DEG and Kj , and Kj 
represent wavevectors at the SBZ corners 1, 2 and 3, 
respectively, in Fig.[T][b). Equivalently, we could express 
the eigenstate as a three-component column vector c = 



(ci,C2,cs) T . Within this basis, the Hamiltonian H, 
to first order in k, is given by H = H + Hi, where 



up 




and 



Hi = hvo k 



H = W 



cos U 
cos (6> k - ^ 




(2) 
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(3) 

Here 9u is the polar angle of the wavevector k from the 
+x direction. 

The eigenvalues of the unperturbed Hamiltonian Hq 
are 



E a = -W, -W, 2W , 



(4) 
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FIG. 4: (color online). Probability densities of the pseudospin 
states in a hexagonal 2DEG superlattice (a) |(r| |)| 2 and (b) 
|(r| J.)| 2 . Note that the amplitudes of the states are localized 
at two different but equivalent sublattices. The centers of the 
potential barrier disks [Fig.fjja)] are shown as 'x' marks, and 
the honeycomb structure is drawn to illustrate the connection 
to the superlattice structure. 



which are also the energies of the states of the superlat- 
tice at k = 0. We now focus on the doubly-degenerate 
eigenstates with eigenvalue —W. We shall find the k- 
dependence of the eigenenergies and eigenvectors of H 
corresponding to these two states within degenerate per- 
turbation theory by treating Hi as a perturbation, which 
is approximate for Hvq k < W. Also, for a wavefunction 
in the form of Eq. {]]) to give a good description of the ac- 
tual wavefunction, the energy to the next planewave state 
[which is h 2 (2K) 2 /2m* - h 2 K 2 /2m*} should be smaller 
than W. Therefore, the approximation is valid within 
the regime Tlvq k < W < ^5—^—, or, equivalently, 



iirh 2 , T „ 8n 2 h 2 
-k <W < 



3m*L 



3m*L 2 ' 



(5) 



where we have used K = 4ir/3L. 

The two eigenvectors of H Q with eigenvalue — W are 



1 




V2 \ _i 



and C2 




(6) 



The term Hi, when restricted to the sub-Hilbert-space 
spanned by the two vectors in Eq. ([6]), is represented by 
a 2 x 2 matrix Hi 



7j— I x " 

2 V 



(7) 



where k x — k cos 9^ and k y = k sin 9^ . After a similarity 
transform M = WHiU with 



we obtain 



1/ l+i -l-i 
2 V -1 + 2 + « 



M = h-£- (k x a x + kytjy 



(8) 



(9) 



Here, a x and a y are the Pauli matrices. Equation ^ is 
just the effective Hamiltonian of graphene^ with a group 
velocity 



v /2 



KK 

2m* 



2nh 
3m* L 



(10) 



Therefore, the group velocity is reduced if the effective 
mass is increased or the lattice parameter of the super- 
lattice is increased, but it is insensitive to the amplitude 
of the external periodic potential, which explains the re- 
sults of the numerical calculations shown in Fig.rjJ Also, 
the size of the linear energy dispersion window [Eq. l[5])] 
is dictated by the value of W, which for the muffin-tin 
potential in Fig. [Ha) with d = 0.663 L is W = 0.172 
in agreement with the numerical calculations^. 
The eigenvalues and the eigenvectors of M are 



E(s, k) = s k , 



and 



kk) 



x/2 (« 



V2' 



(11) 



(12) 



respectively, where s = ±1 is a band indes^. The vectors 



are the up and the down pseudospin 



o ) and ( 1 , 

eigenstates of a z , respectively^. 

The up and the down pseudospin eigenstates may be 
expressed within the basis of the original Hamiltonian H 
using Eqs. (jHJ) and © as 



and 
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II) 
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respectively 2 -. The real space pseudospin wavefunctions 
(r| t) and (r| J.) are obtained by putting the coefficients 
in Eqs. {32]) and dHJ) into Eq. Figures (Ha) andS^b) 
show |(r| |)| 2 and |(r| |)| 2 , respectively. Note that the 
up and the down pseudospin states are seen as localized 
at one of the two equivalent sublattices formed by the 
external periodic potential, in perfect analogy with the 
behavior in graphene. 

Let us now consider the Landau levels for the above 
hexagonal 2DEG superlattice in a magnetic field B = B z 
when the Fermi level is at the Dirac point energy. In ex- 
act analogy with graphene, the low-energy Landau levels 
shifted by the energy W (i. e. , having the energy zero at 
the Dirac point energy) are E n = huj c sgn(n) yf\n~\, where 
uj c = vq y/\e\B/2hc is the cyclotron frequency^^. Here, 
for an appropriately constructed superlattice potential, 
the half-integer quantum Hall effect^ should be observ- 
able. 

In conclusion, we have shown that chiral massless Dirac 
fcrmions are generated if an appropriate nanometer-scale 
periodic potential with hexagonal symmetry is applied 
to a conventional 2DEG in semiconductors. These quasi- 
particles have a linear energy dispersion, with a group 
velocity half the value of the states before applying the 
periodic potential, and a wavefunction whose chiral struc- 
ture exactly the same as that of graphene. The up and 



the down pseudospin states are shown to be localized 
at two different but equivalent sublattices formed by the 
superlattice potential. The quasiparticle group velocity 
moreover is tunable by changing the effective mass of 
the original 2DEG or the lattice parameter of the super- 
lattice potential. Our findings thus provide a new class 
of systems for experimental investigations and practical 
applications of 2D massless Dirac quasiparticles. 
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Note added in proof. - After we finished writing up 
this manuscript, we became aware of Ref4^ in which a 
similar methodology to the one developed independently 
here was developed. The objects of investigation (which 
are cold fermionic atoms confined in honeycomb optical 
lattices) and the main ideas and conclusions in Ref. 31 , 
however, are qualitatively different from those in this Let- 
ter. 



* Electronic address: sglouie@berkeley.edu 

1 K. S. Novoselov, D. Jiang, F. Schedin, T. Booth, V. V. 
Khotkevich, S. V. Morozov, and A. K. Geim, Proc. Natl. 
Acad. Sci. USA. 102, 10451 (2005). 

2 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
F. A. A., Nature 438, 197 (2005). 

3 Y. Zhang, J. W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

4 C. Berger, Z. Song, X. Li, X. Wu, N. Brown, P. N. First, 
and W. A. de Heer, Science 312, 1191 (2006). 

5 P. R. Wallace, Phys. Rev. 71, 622 (1947). 

6 T. Ando and T. Nakanishi, J. Phys. Soc. Jpn. 67, 1704 
(1998). 

7 T. Ando, F. Nakanishi, and R. Saito, J. Phys. Soc. Jpn. 
67, 2857 (1998). 

8 J. W. McClure, Phys. Rev. 104, 666 (1956). 

9 M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. 
Phys. 2, 620 (2006). 

10 P. L. McEuen, M. Bockrath, D. H. Cobden, Y.-G. Yoon, 
and S. G. Louie, Phys. Rev. Lett. 83, 5098 (1999). 

11 C.-H. Park, L. Yang, Y. W. Son, M. L. Cohen, and S. G. 
Louie, Nature Phys. 4, 213 (2008). 

12 C.-H. Park, Y.-W. Son, L. Yang, M. L. Cohen, and S. G. 
Louie, Nano Lett. 8, 2920 (2008). 

13 C.-H. Park, L. Yang, Y.-W. Son, M. L. Cohen, and S. G. 
Louie, Phys. Rev. Lett. 101, 126804 (2008). 

14 G. Grynberg, B. Lounis, P. Verkerk, J.-Y. Courtois, and 
C. Salomon, Phys. Rev. Lett. 70, 2249 (1993). 

15 S.-L. Zhu, B. Wang, and L.-M. Duan, Phys. Rev. Lett. 98, 



260402 (2007). 

16 C. Wu, D. Bergman, L. Balents, and S. D. Sarma, Phys. 
Rev. Lett. 99, 070401 (2007). 

17 C. Wu and S. D. Sarma, Phys. Rev. B 77, 235107 (2008). 

18 L. H. Haddad and L. D. Carr, arXiv:0803.3039vl. 

19 L. B. Shao, S.-L. Zhu, L. Sheng, D. Y. Xing, and Z. D. 
Wang, Phys. Rev. Lett. 101, 246810 (2008). 

20 J. F. Devreese and F. M. Peeters, eds., The Physics of 
the Two-Dimensional Electron Gas (Plenum, New York, 
1987). 

21 E. Ribeiro, E. Miiller, F. Heinzel, H. Auderset, K. Ensslin, 
G. Medeiros-Ribeiro, and P. M. Petroff, Phys. Rev. B 58, 
1506 (1998). 

22 E. Ribeiro, R. D. Jaggi, T. Heinzel, K. Ensslin, 
G. Medeiros-Ribeiro, and P. M. Petroff, Phys. Rev. Lett. 
82, 996 (1999). 

23 E. Ribeiro, R. Jaggi, T. Heinzel, K. Ensslin, F. G. 
Medeiros-Ribeiro, and P. M. Petroff, Microelectron. Eng. 
61, 674 (1999). 

24 P. Y. Yu and M. Cardona, Fundamentals of Semiconduc- 
tors: Physics and Materials Properties (Springer, New 
York, 2001). 

25 The Fourier component W, defined to be the value of the 
external periodic potential [Fig. QJa)] in Fourier space at 
the smallest reciprocal lattice vector, is given by W = 

Y^J\ ^~7fjM Uo where Ji(x) is the first-order Bessel func- 
tion. For a given Uo, the maximum value of W is W = 
0.172 Uo at d — 0.663 L. In our numerical calculations, 
we take advantage of this. For comparison, if d — 0.5 L, 
W = 0.145 U - 



5 



C. Albrecht, J. H. Smet, D. Weiss, K. von Klitzing, R. Hen- 
nig, M. Langenbuch, M. Suhrke, U. Rossler, V. Umansky, 
and H. Schweizer, Phys. Rev. Lett. 83, 2234 (1999). 
In fabricated 2DEG superlattices, small fluctuations will 
inevitably occur in the amplitude and periodicity of 
the potential. These have been studied theoretically 
(Refs . 32 ' 33 ' 34 ) for systems with honeycomb symmetry (like 
graphene) and found to be of little effect on the DOS and 
transport properties at low concentrations. In particular, a 
recent ab initio study (Ref.— ) of graphene shows insignifi- 
cant changes in the linear behavior in the DOS around the 
Dirac point energy up to 5 % in concentration of boron 
and nitrogen impurities. We expect the same robustness 
should occur for the 2DEG superlattices. 
Similar procedures can be applied to find the energy and 
the wavefunction of massless Dirac fermions with wavevec- 
tors near K' [Fig. \Qb)]. 
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a phase factor, the eigenvector in the original planewave 
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sin(% + f) 



( S = -l) 



The expression for the (relativistic) Landau levels is valid 
when the spacing between Landau levels are much smaller 
than the size of the energy window within which the disper- 
sion is linear, regardless of the existence of higher energy 
states, just as in graphene. For example, in a triangular 
2DEG superlattice with m* = 0.02 m e , U = 200 meV, 
L = 20 nm and d = 13.3 nm, this condition is equivalent 
to B < 2 T. 
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